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Abstract: Wc define a new theory of discrete Riemann surfaces and present 
its basic results. The key idea is to consider not only a cellular decomposition 
of a surface, but the union with its dual. Discrete holomorphy is defined by a 
straightforward discretisation of the Cauchy- Riemann equation. A lot of classical 
results in Riemann theory have a discrete counterpart, Hodge star, harmonicity, 
Hodge theorem, Weyl's lemma, Cauchy integral formula, existence of holomor- 
phic forms with prescribed holonomies. Giving a geometrical meaning to the 
construction on a Riemann surface, we define a notion of criticality on which 
we prove a continuous limit theorem. We investigate its connection with crit- 
icality in the Ising model. We set up a Dirac equation on a discrete universal 
spin structure and we prove that the existence of a Dirac spinor is equivalent to 
criticality. 
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1. Introduction 

We present here a new theory of discrete analytic functions, generalising to 
discrete Riemann surfaces the notion introduced by Lelong-Ferrand |LFj . 

Although the theory defined here may be applied wherever the usual Rie- 
mann Surfaces theory can, it was primarily designed with statistical mechanics, 
and particularly the Ising model, in mind [McC WiliD] . Most of the results can 
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be understood without any prior knowledge in statistical mechanics. The other 
obvious fields of application in two dimensions are electrical networks, elasticity 
theory, thermodynamics and hydrodynamics, all fields in which continuous Rie- 
mann surfaces theory gives wonderful results. The relationship between the Ising 
model and holomorphy is almost as old as the theory itself. The key connection 
to the Dirac equation goes back to the work of Kaufman [K] and the results in 
this paper should come as no surprise for workers in statistical mechanics; they 
knew or suspected them for a long time, in one form or another. The aim of 
this paper is therefore, from the statistical mechanics point of view, to define a 
general theory as close to the continuous theory as possible, in which claims as 
"the Ising model near criticality converges to a theory of Dirac spinors" are given 
a precise meaning and a proof, keeping in mind that such meanings and proofs 
already exist elsewhere in other forms. The main new result in this context is 
that there exists a discrete Dirac spinor near criticality in the finite size Ising 
model, before the thermodynamic limit is taken. Self-duality, which enabled 
the first evaluations of the critical temperature [KWilOnsllWanSO] . is equivalent 
to criticality at finite size. It is given a meaning in terms of compatibility with 
holomorphy. 

The first idea in order to discretise surfaces is to consider cellular decompo- 
sitions. Equipping a cellular decomposition of a surface with a discrete metric, 
that is giving each edge a length, is sufficient if one only wants to do discrete har- 
monic analysis. However it is not enough if one wants to define discrete analytic 
geometry. The basic idea of this paper is to consider not just the cellular decom- 
position but rather what we call its double, i.e. the pair consisting of the cellular 
decomposition together with its Poincare dual. A discrete conformal structure is 
then a class of metrics on the double where we retain only the ratio of the lengths 
of dual edgesQ. In Ising model terms, a discrete conformal structure is nothing 
more than a set of interaction constants on each edge separating neighbouring 
spins in an Ising model of a given topology. 

A function of the vertices of the double is said to be discrete holomorphic if 
it satisfies the discrete Cauchy-Riemann equation, on two dual edges {x,x') and 

{y,y'), 

f{y')-f{y) ^ J{^')-f{^) 

£iy,y') ' eix,x') ' 

This definition gives rise to a theory which is analogous to the classical theory 
of Riemann surfaces. We define discrete differential forms on the double, a Hodge 
star operator, discrete holomorphic forms, and prove analogues of the Hodge 
decomposition and Weyl's lemma. We extend to our situation the notion of pole 
of order one and we prove existence theorems for meromorphic differentials with 
prescribed poles and holonomies. Similarly, we define a Green potential and a 
Caucliy integral formula. 

Up to this point, the theory is purely combinatorial. In order to relate the 
discrete and continuous theories on a Riemann surface, we need to impose an 
extra condition on the discrete conformal structure to give its parameters a geo- 
metrical meaning. We call this semi-criticality in Sect.[3l The main result here is 



^ By definition, a discrete Riemann surface is a discrete surface equipped with a discrete 
conformal structure in this sense. 
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Fig. 1: The discrete Cauchy-Riemann equation. 



that the hmit of a pointwise convergent sequence of discrete holomorphic func- 
tions, on a refining sequence of semi-critical cellular decompositions of the same 
Riemann surface, is a genuine holomorphic function on the Riemann surface. If 
one imposes the stronger condition of criticality on the discrete conformal struc- 
ture, one can define a wedge product between functions and 1-forms which is 
compatible with holomorphy. 

Finally, for applications of this theory to statistical physics, one needs to 
define a discrete analogue of spinor fields on Riemann surfaces. In Sect. |4] we 
first define the notion of a discrete spin structure on a discrete surface. It sheds 
an interesting light onto the continuous notion, allowing us to redefine it in 
explicit geometrical terms. In the case of a discrete Riemann surface we then 
define a discrete Dirac equation, generalising an equation appearing in the Ising 
model, and show that criticality of the discrete conformal structure is equivalent 
to the existence of a local massless Dirac spinor field. 

In Sect. [2 we present definitions and properties of the theory which are purely 
combinatorial. First, in the empty boundary case, we recall the definitions of dual 
cellular complexes, notions of dcRham cohomology. We define the double A, we 
present the discrete Cauchy-Riemann equation, the discrete Hodge star on A, the 
Laplacian and the Hodge decomposition. In Subsect. l2.2[ we prove Dirichlet and 
Neumann theorems, the basic tools of discrete harmonic analysis. In Subsect. 12.31 
we prove existence theorems for 1-forms with prescribed poles and holonomies. 
In Subsect. 12.41 we deal with the basic difficulty of the theory: The Hodge star 
is defined on A while the wedge product is on another complex, the diamond 
(}, obtained from F or F* by the procedure of tile centering |GS87| . We prove 
Weyl's lemma. Green's identity and Cauchy integral formulae. 

In Sect. [3l we define semi-criticality and criticality and prove that it agrees 
with the usual notion for the Ising model on the square and triangular lattices. 
We present Voronoi and Delaunay semi-critical maps in order to give examples 
and we prove the continuous limit theorem. We prove that every Riemann surface 
admits a critical map and give examples. On a critical map, the product between 
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functions and 1-forms is compatible with holomorphy and yields a polynomial 
ring, integration and derivation of functions. We give an example showing where 
the problems are. 

In Sect.m we set up the Dirac equation on discrete spin structures. We moti- 
vate the discrete universal spin structure by first showing the same construction 
in the continuous case. We show discrete holomorphy property for Dirac spinors, 
we prove that criticality is equivalent to the existence of local Dirac spinors and 
present a continuous limit theorem for Dirac spinors. 

2. Discrete Harmonic and Holomorphic Functions 

In this section, we are interested in properties of combinatorial geometry. The 
constructions are considered up to homeomorphisms, that is to say on a combi- 
natorial surface, as opposed to Sect. [3]whcre criticality implies that the discrete 
geometry is embedded in a genuine Riemann surface. 

2.1. First definitions. Let S be an oriented surface without boundary. A cellular 
decomposition of is a partition of S into disjoint connected sets, called cells, 
of three types: a discrete set of points, the vertices Iq; a set of non intersecting 
paths between vertices, the edges Fi; and a set of topological discs bounded by 
a finite number of edges and vertices, the oriented faces F2. A paramctrisation 
of each cell is chosen, faces are mapped to standard polygons of the euclidean 
plane, and edges to the segment (0, 1); we recall particularly that for each edge, 
one of its two possible orientations is chosen arbitrarily. We consider only locally 
finite decompositions, i.e. any compact set intersects a finite number of cells. In 
each dimension, we define the space of chains Ct{F) as the Z-module generated 
by the cells. The boundary operator d : Ck{F) —f Ck-i{F) partially encodes the 
incidence relations between cells. It fulfills the boundary condition dd = 0. 

Wc now describe the dual cellular decomposition F* of a cellular de- 
composition r' of a surface without boundary. We refer to |Veb| for the general 
definition. Though we formally use the paramctrisation of each cell for the def- 
inition of the dual, its combinatorics is intrinsically well defined. To each face 
F G F2 we define the vertex F* £ Fq inside the face F, the preimage of the ori- 
gin of the euclidean plane by the paramctrisation of the face. Each edge e £ Fi, 
separates two faces, say Fi,F2 G A (which may coincide), hence is identified 
with a segment on the boundary of the standard polygon corresponding to Fi, 
respectively F2 . We define the dual edge e* <E F^ as the preimage of the two seg- 
ments in these polygons, joining the origin to the point of the boundary mapped 
to the middle of e. It is a simple path lying in the faces Fi and -F2, drawn be- 
tween the two vertices F* and F2 (which may coincide), cutting no edge but e, 
once and transversely. As the surface is oriented, to the oriented edge e we can 
associate the oriented dual edge e* such that (e, e*) is direct at their crossing 
point. To each vertex v £ /q, with vi, . . . ,Vn € Fq as neighbours, we define the 
face V* <E F2 by its boundary dv* = {v, vi)* + . .. + {v, Vk)* + . . . + {v, w„)*. 

Remark L F* is a. cellular decomposition of U |Veb| . If we choose a paramctrisa- 
tion of the cells of F* , we can consider its dual F**; it is a cellular decomposition 
homeomorphic to F but the orientation of the edges is reversed. The bidual of 
e S A is the reversed edge e** = — e (see Fig. [2|). 
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0. The vertex dual to a face. 1. Dual edges. 2. The face dual to a vertex. 



Fig. 2: DuaUty. 

The double ^4 of a cellular decomposition is the union of these two dual 
cellular decompositions. We will speak of a fc-cell of yl as a k-ce\l of either F or 

A discrete metric i is an assignment of a positive number £{e) to each edge 
e £ Ai, its length. For convenience the edge with reversed orientation, — e, will 
be assigned the same length: £{—e) := £{e). Two metrics £,£' : A\ —* (0, +oo) 
belong to the same discrete conformal structure if the ratio of the lengths 
p[e) := = 77^, on each pair of dual edges e S A, e* € are equal. 

A function / on yl is a function defined on the vertices of F and of F* . Such 
a function is said to be holomorphic if, for every pair of dual edges {x,x') G Fi 
and {y,y') = {x,x')* e F^ , it fulfills 

fiy')-.fiy) _ J{^')-f{^) 

£iy,y') ' £{x,x') 

It is the naive discretisation of the Cauchy-Riemann equation for a function /, 
which is, in local orthonormal coordinates {x,y): 

dy dx 

Here, we understand two dual edges as being orthogonal. 

This equation, though simple, was never considered in such a generality. It 
was introduced by Lelong-Ferrand [LFj for the decomposition of the plane by the 
standard square lattice 1? . It is also called monodiffric functions; for background 
on this topic, see |Duf] . Polynomials of degree two, restricted to the square lattice, 
give examples of monodiffric functions. See also the works of Kenyon [Ken] and 
Schramm and Benjamini [BS96| who considered more than lattices. 

The usual notions of deRham cohomology arc useful in this setup. We said 
that /c-chains are elements of the Z- module Ck{A), generated by the fc-cells, its 
dual space C'^(yl) := Hom (Cfc(yl),C) is the space of fc-cochains. We will denote 
the coupling by the usual integral and functional notations: /(x) for a function 
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/ e C"(yl) on a vertex x G ylo; a for a 1-form a G C^(yl) on an edge e G Ai; 
and JJpOj for a 2-forni a; e (A) on a face F g . 

The coboundary d : C'^(yl) C''~^^{A) is defined by the Stokes formula 
(with the same notations as before): 



df / x')) = fix') - fix) jj da:= (ta. 

(x.x') F OF 

As the boundary operator splits onto the two dual complexes F and -T*, the 
coboundary d also respects the direct sum C^iA) ~ C'^(r) © C^iF*). 

The Cauchy-Riemann equation can be written in the usual form *df ~ ^idf 
for the following Hodge star * : C''iA) C^~''(yl) defined by: 



a —pie*) / a 

J e* 

We extend it to functions and 2-forms by: 

*f := fiF*), M^) := 



As, by definition, for each edge e & Ai, p(e)p(e*) = ^(^j^ = 1, 

the Hodge star fulfills on fc-forms, *^ — (— l)'^(^^'^'ldcfc(/i)- 
It decomposes 1-forms into —i, respectively +i, eigenspaces, called type 
(1,0), rcsp. type (0,1) forms: 

C\A)^C^^'°\A)(BC<'^'^HA). 
The associated projections arc denoted: 

7r(i,o) -i(Id + z*):Ci(^)^C(i'0)(yl), 

7r^o.i) = \ild-t*):C\A)^C^'''\A). 

A 1-form is holomorphic if it is closed and of type (1,0): 

a G i7^(yl) da ~ and * a = —ia. 

It is meromorphic with a pole at a vertex a; G ylo if it is of type (1,0) and not 
closed on the face x*. Its residue at x is defined by 



ReSa;(Q;) := * 



The residue theorem is merely a tautology in this context. 

We define d', d" , the composition of the coboundary with the projections on 
eigenspaces of * as its holomorphic and anti-holomorphic parts: 

r(i,o) ° d, d" := 7r(o.i) o d 
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d' :— do TT^Q 2) , d" :— do tt^^ 



0) 



from 1-forms to 2-forms and d' = rf" = on 2-forms. They verify d'^ = and 
d"^ = 0. 

The usual discrete laplacian, which sphts onto F and 7^* independently, 
reads A := —d *d* — *d*da.s expected. Its formula for a function / G C'^{A) 
on a vertex x G Aq, with xi, . . . ,Xn as neighbours is 

n 

{Af){x) = J2 '^(^' (/(^) - /(^fc)) ■ (2.1) 
fe=i 

As in the continuous case, it can be written in terms of d' and d" operators: For 
functions, A = i* [d'd" — d"d'), in particular holomorphic and anti-holomorphic 
functions are harmonic. The same result holds for 1-forms. 

In the compact case, the operator d* = — * is the adjoint of the coboundary 
with respect to the usual scalar product, (/,(?) :— J^xgAq /(^)ff(^) '^^ functions, 
similarly on 2-forms and 



(3 I on 1— forms. 



It gives rise to the Hodge decomposition. 

Proposition 1 (Hodge theorem). In the compact case, the k-forms are de- 
composed into orthogonal direct sums of exact, coexact and harmonic forms: 



C'^iA) = Im d Im d* ®^ Ker A, 



harmonic forms are the closed and coclosed ones: 



Ker A = Ker d n Ker d* . 



In particular the only harmonic functions are locally constant. Harmonic 1-forms 
are also the sum of holomorphic and anti-holomorphic ones: 



Ker A ^ Ker d' 0^ Ker d" . 



Beware that A being disconnected, the space of locally constant functions is 
2-dimensional. The function e which is +1 on F and —1 on F* is chosen as the 
second basis vector. 

The proof is algebraic and the same as in the continuous case. As the Lapla- 
cian decomposes onto the two dual graphs, this result tells also that for any 
harmonic 1-form on F, there exists a unique harmonic 1-form on the dual graph 
F* such that the couple is a holomorphic 1-form on A, it's simply a^. := i* ar- 
These decompositions don't hold in the non-compact case; there exist non-closed 
and/or non-co-closed, harmonic 1-forms. 
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2.2. Dirichlet and Neumann problems. 

Proposition 2 (Dirichlet problem). Consider a finite connected graph F, 
equipped with a function p on the edges, and a certain non-empty set of points 
D marked as its boundary. For any boundary function f^ : {dF)o — > C, there 
exists a unique function f , harmonic on Fq\ D such that f\dr — f^- 

Wc refer to the usual laplacian defined by Eq. 12. f I 

If z= 0, the solution is the null function. Otherwise, it is the minimum of 
the strictly convex, positive functional / {df,df), proper on the non-empty 
affine subspace of functions which agree with f^ on the boundary. □ 

Definition 1. Given F a cellular decomposition of a compact surface with bound- 
ary S, define the double := SU S, union with the opposite oriented surface, 
along their boundary. The double F^ is a cellular decomposition of the compact 
surface E^. Consider its dual F^* and define F* := SCiF'^* . We don't take into 
account the faces of F^* which are not completely inside S but we do consider 
the half-edges dual to boundary edges of F as genuine edges noted (dF*)i and 
define {dF*)o := 7^;^^* n dS as the set of their boundary vertices. 

A function p on the edges of F yields an extension to /\* by defining p{e*) :— 

pR- 

Remark 2. F* is not a cellular decomposition of the surface; the half-edges dual 
to boundary edges do not bound any face of F*. 

Proposition 3 (Neumann problem). Consider F a cellular decomposition 
of a disk, equipped with a function p on its edges. Choose a boundary vertex 
yo G {dF*)Q, a value /o G C, and on the set of boundary edges e € {dF*)i, not 
incident to yo, a 1-form a. 

Then there exists a unique function f, harmonic on F* \ {dF*)^ such that 
fivo) = /o CLnd J^df = ^^a for all e G {dF*)i not incident to yo- 

It is a dual problem. Let eg S {dF*)i, be the edge incident to yo and cq G 
{dF)i its dual. Consider, on the set of boundary edges e G {dF)i different from 
eg, the 1-form defined by i * a. Integrating it along the boundary, we get a 
function on {dF)o, well defined up to an additive constant. The Dirichlet 
theorem gives us a function / harmonic on Iq \ {dF)o corresponding to 
Integrating the closed 1-form i ^ df on F* yields the desired harmonic function 
/•□ 

Remark 3. The number of boundary points in F is the same as in F* , and as every 
harmonic function on _r, when the surface is a disk, defines a harmonic function 
on F* such that their couple is holomorphic, unique up to an additive constant, 
the space of holomorphic functions, resp. 1-forms, on the double decomposition 
with boundary A is |(9yl)o|/2 -I- 1, resp. |(9yl)o|/2 — 1 dimensional. 

The theorem is true for more general surfaces than a disk but the proof is 
different, see the author's PhD thesis fM^. There are i"^ versions of these theorems 
too. 
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2.3. Existence theorems. We have very similar existence theorems to the ones in 
the continuous case. We begin with the main difference: 

Proposition 4. The space oj discrete holomorphic 1-forms on a compact surface 
without boundary is of dimension twice the genus. 

The Hodge theorem imphes an isomorphism between the space of harmonic 
forms and the cohomology group of A. It is the direct sum of the cohomology 
groups of r and of F* and each is isomorphic to the cohomology group of the 
surface which is 2g dimensional on a genus g surface. It splits in two isomorphic 
parts under the type (1,0) and type (0,1) sum. As any holomorphic form is 
harmonic, the dimension of the space of holomorphic 1-forms is 2g.O 

We can give explicit basis to this vector space as in the continuous case |Sie| . 
To construct them, we begin with meromorphic forms: 

Proposition 5. Let S be a compact surface with boundary. For each vertex 
X Ao\ 9<C^, and a simple path X on A going from x to the boundary there exists 
a pair of meromorphic 1-forms ax^Px with a single pole at x, with residue +1 
and which have pure imaginary, respectively real holonomies, along loops which 
don't have any edge dual to an edge of X. 

Proposition 6. Let S be a compact surface. For each pair of vertices x, x' G Aq 
with a simple path X on A from x to x' , there exists a unique pair of meromor- 
phic 1-forms Oix,x' , Px,x' with only poles at x and x' , with residue +1 and — 1 
respectively, and which have pure imaginary, respectively real holonomies, along 
loops which don't have any edge dual to an edge of X. 

In both cases, the forms are (Id + i*)df with / a solution of a Dirichlet 
problem at x (and x') for a and of a Neumann problem on the surface split open 
along the path A for (3. The uniqueness in the second proposition is given by 
the difference: the poles cancel out and yield a holomorphic 1-form with pure 
imaginary, resp. real holonomies, so its real part, resp. imaginary part, can be 
integrated into a harmonic, hence constant function. So this part is in fact null. 
Being a holomorphic 1-form, the other part is null too. We refer to the author's 
PhD thesis [M] for details.D 

As in the continuous case, it allows us to construct holomorphic forms with 
(no poles and) prescribed holonomies: 

Corollary 1 Let A,B ^ Zi{A) be two non-intersecting simple loops such that 
there exists exactly one edge of A dual to an edge of B (dual loops). There exists 
a unique holomorphic 1-form 'Pab such that lic{Jg<pAB) = 1 ^'^'^ J^'^AB G 
for every loop 7 which doesn't have any edge dual to an edge of A. 

We decompose A in two paths and Xy. It gives us two 1-forms f3x,y and 
(3y^x, then 

'^AB-= l^{f3.,y+Py,.) (2.2) 
ZlTT 

fulfills the conditions. □ 
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2.4- The diamond <0 and its wedge product. Following |Whit| . we define a wedge 
product, on another complex, the diamond (}, constructed out of the double A: 
Each pair of dual edges, say {x,x') G A and {y,y') = {x,x')* G F*, defines (up 
to homeomorphisms) a four-sided polygon {x,y, x' ,y') and all these constitute 
the faces of a cellular complex called (see Fig. [3]). 



.■9.. y' 



X 



X 



Fig. 3: The diamond (}. 



On the other hand, from any cellular decomposition of a surface by four- 
sided polygons one can reconstruct the double A. A difference is that A may 
not be disconnected in two dual pieces F and _r*, it is so if each loop in is 
of even length; we will restrict ourselves to this simpler case. This is not very 
restrictive because from a connected double, refining each quadrilateral in four 
smaller quadrilaterals, one gets a double disconnected in two dual pieces. 



Definition 2. A discrete surface with boundary is defined by a quadrilateral 
cellular decomposition <0> of an oriented surface with boundary such that its double 
complex A is disconnected in two dual parts. 



This definition is a generalisation of the more natural previous Definition [TJ 
It allows us to consider any subset of faces of ^ as a domain yielding a discrete 
surface with boundary. While any edge of A has a dual edge, a vertex of A 
has a dual face if and only if it is an inner vertex. Punctured surfaces can be 
understood in these terms too: An inner vertex v € Aq is a puncture if it is 
declared as being on the boundary and its dual face v* removed from A2. 

We construct a discrete wedge product, but while the Hodge star lives on the 
double A, the wedge product is defined on the diamond {): A : C'^{(^) x C'({)) 
C'=+'(0). It is defined by the following formulae, for f,ge C"{<}), a,/3 G C\<}) 
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if ■ 9){x) ■■=f{x) ■ g{x) forxe^o, 

f.a:=M±I^ f {x,y)a for (.t, e Oi, 

-"-it / °/"- / 

{xi,X2,X3,X4) ' {Xk-l,Xk){Xk,Xk + l) {x k + 1 ,X k) {x k ,X ^ - l) 

f{xi)+f{x2)+f{x3)+f{Xi) f f 

/..:= JJ. 

{xi,X2,X3,Xi) 3:2, 2:3, 3:4) 

for (xi,a;2,X3,a;4) G 02- 

Lemma 1. is a derivation with respect to this wedge product. 

To take advantage of this property, one has to relate forms on (} and forms 
on A where the Hodge star is defined. Wc construct an averaging map A 
from C*(<^) to C*{A). The map is the identity for functions and defined by the 
foUowing formulae for 1 and 2-forms: 



,x') \x,y) {y.x'){x,y') {y',x')/ 



(2.3) 



^' ^ ^(xk-Vk,x,yk-i) 

where notations are made clear in Fig. |4l With this definition, d/iA = Ad^, but 
the map A is neither injcctive nor always surjective, so we can neither define a 
Hodge star on (} nor a wedge product on A. An element of the kernel of A is 
given for example by d^^e, where £ is +1 on F and —1 on 7^*. 

On the double A itself, we have pointwise multiplication between functions, 
functions and 2-forms, and we construct an heterogeneous wedge product for 
1-forms: with a, /3 G C^{A), define a A /3 G Ci«» by 



a/\P := J a J (3 + J a J p. 

{x,y,x',y'} (x,x')(y,y') (y,y')(x',x) 

It verifies A(a^) A A{I3^) — A /3^, the first wedge product being between 
1-forms on A and the second between forms on <). Of course, we also have for 
integrable 2-forms: 

" // "^^^"^^ ^11 ^^"^^^ ^^11 
O2 ^2 A2 

And for a function /, 

A2 A2 
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(EJD (El 

Fig. 4: Notations. 



whenever / • lo^ is integrable. 

Explicit calculation shows that for a function / G (7°(yl), denoting by Xx the 
characteristic function of a vertex .t e tIq, = ~//^^ / • *^XrE- So by 

linearity one gets Weyl's lemma: a function / is harmonic iff for any compactly 
supported function g G C^{A), 



II 



One checks also that the usual scalar product on compactly supported forms 
on A reads as expected: 



a A */3. 



In some cases, for example, the decomposition of the plane by lattices, the av- 
eraging map A is surjective. We define the inverse map B : [A) C^{^) /Ker A 
and Z\<) := d^B ^ d and we then have 

Proposition 7 (Green's identity). For two functions f,g on a compact do- 
main D d ()>2- 

if-A^g-g- A^f) - i [f ■ B*dg - g ■ B*df) = 0. 

D JdD 

This means that for any representatives of the classes in C"'^(<0>)/Ker A the 
equality holds, but each integral separately is not well defined on the classes. 
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2. 5. Cauchy integral formula. 

Proposition 8. Let A a double map and D a compact region of <^2 homeomor- 
phic to a disc. Consider an interior edge {x, y) € D; there exists a meromorphic 
1-form i'x,y G C'^{D\ {x,y)) such that the holonomy §^Vx.y along a cycle 7 in 
D only depends on its homology class in D\ (x, y), and v^^y ~ 2iiT . 

Consider the meromorphic 1-form pLx.y ^ + cty e C^{An D) defined by 
existence Theorcm[5]on D. It is uniquely defined up to a global holomorphic form 
on D. Its only poles are x and y of residue +1 so it verifies a similar holonomy 
property, but on D\{x,y). We define a 1-form v^.y on <^ n D \ i?, such that 
l^x,y = Avx.y in the following way: Let J^^ v^^y ■= A, a fixed value, and for an 
edge ix',y') £ Di, with x' e Tq, y' e Tq*, given two paths in £>, \%, € Ci(r) 
and \y e Ci(r*) respectively from x' to x and from y to y', 

/ ■ — / l^x^y ~t- / ^a:,y ~t~ / l^x.y y l^x,y-! 

J{x',y') JXI, J(x.A) Jxl J[7] 

where [7] is the homology class of X^, -f (x, y) + -t- (y', x') on the punctured 
domain. □ 

i'x,y is the discrete analogue of -^r^ with zq = {x, y). It is closed on every face 
of D\R. By definition, the average of i>x,y on the double map is the meromorphic 

form Aiyx,y = tix.y 

It allows us to state 

Proposition 9 (Cauchy integral formula). Let D be a compact connected 

.subset of <^2 and {x, y) e Di two interior neighbours of D with a non-empty 
boundary. For each function f G C*'(yl), 

/ • J^x,y = / / a f A flx.y + 2i7r . 

dD J J D ^ 

The proof is straightforward: The edge (a;,y) boimds two faces in £), let 
R = (abed) the rectangle made of these faces (see Fig. [5]). 




, y b 




Fig. 5: The rectangle i? in a domain D defined by an edge (x, y) E (^i. 



On D\R, 

rfo(/ • '^x,y) = d^f A Ux.y + f ■ d^Vx,y 
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The (1,0) part of df disappears in the wedge product against the holomorphic 
form fJ-x.y, so we can substitute 

rfo/ A Uj;^y = A Aiy^^y = d"f A ^i.j,^y. 

Integrating over D, closed on D \ R, we get: 



dD J J D\R JdR 



Exphcit calculus shows that §qj^ f ■ Ux,y ~ jj^j d" f A ^x,y + 2in ^^^''^^^^' . □ 

Remark 4- Since for all a £ C^(<)), the locally constant function e defined by 
e{r) = +1, e{r*) = — 1, verifies £ • a = 0, an integral formula will give the same 
result for a function / and / + Ae. Therefore such a formula can not give access 
to the value of the function at one point but only to its average value at an edge 
oiO. 

Corollary 2 For f g f2{ A) a holomorphic function, the Cauchy integral formula 
reads, with the same notations, 

fix) + fiy) _ 1_ 

IdD 



2 ^2^^'P 



The Green function on the lattices (rectangular, triangular, hexagonal, Kagome, 
square/octogon) is exactly known in terms of hyperelliptic functions ( |Hug| and 
references in Appendix 3). As the potential is real, it means that the discrete 
Dirichlet problem on these lattices can be exactly solved this way, once the 
boundary values on the graph and its dual are given: if these values are real 
and r and imaginary on its dual, the solution is real on F and pure imaginary 
on the dual so the contributions f{x) and f{y) are simply the real and imag- 
inary parts of the contour summation respectively. Unfortunately, this pair of 
boundary values arc not independant but related by a Dirichlet to Neumann 
problem |CdV96j . 



3. Criticality 

The term criticality, as well as our motivation to investigate discrete holomorphic 
functions, comes from statistical mechanics, namely the Ising model. A critical 
temperature is defined that restrains the interaction constants, interpreted here 
as lengths. We will see these geometrical constraints in Sect. 13.31 

Technically, as far as the continuous limit theorem is concerned, a weaker 
property, called semi- criticality is sufficient, it gives us a product between func- 
tions and forms. Moreover, at criticality, this product will be compatible with 
holomorphy. 
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3.1. Semi-cnticality. Define Ce := {(r,t) : r > 0,< e R/6'Z}/(0,i) ~ (0,f) with 
the metric ds^ := dr^ + r'^dt^ as the standard cone of angle 9 > [Tro| . 

The cones can be realized by cutting and pasting paper, demonstrating their 
local isometry with the euclidean complex plane. 

Let S he a, compact Riemann surface and P <Z S & discrete set of points. A 
flat Riemannian metric with P as conic singularities is an atlas {Zij^ : 
Ux ^ U' C Ce^ >}x£P of open sets Ux C S, a, neighbourhood of a singularity 
X G P, into open sets of a standard cone, such that the singularity is mapped 
to the vertex of the cone and the changes of coordinates Cu,v : U OV C are 
euclidean isometrics. 

There is a lot of freedom allowed in the choice of a flat metric for a given 
closed Riemann surface S: Any finite set P of points on X! with a set of angles 
9x > for every x G P such that 2Trx{S) = X^^gp^^^ ~ ^^)' defines uniquely a 
Riemannian flat metric on E with these conic singularities and angles [Troj . 

Consider such a flat riemannian metric on a compact Riemann surface S and 
{A, i) a double cellular decomposition of S as before. 

Definition 3. {A,£) is a semi-critical map for this flat metric if the conic 
singularities are among the vertices of A and <^ can be realized such that each face 
[x, y, x' , y') € <0>2 is mapped, by a local isometry Z preserving the orientation, to a 
four-sided polygon {Z{x), Z{y), Z{x'), Z(y')) of the euclidean plane, the segments 
[Z{x), Z{x')] and [Z{y), Z{y')] being of lengths l{x,x'), £{y,y') respectively and 
forming a direct orthogonal basis. We name S{A,£) the supremum of the lengths 
of the edges of <C> . 

The local isometric maps Z are discrete holomorphic. 

Voronoi and Delaunay complexes [PS85| are interesting examples of 
semi-critical dual complexes. Any discrete set of points Q on a flat Riemannian 
surface, containing the conic singularities, deflnes such a pair: 

We first define two partitions V and D oi E into sets of three types: 2-sets, 
1-sets and 0-sets, and then show that they are in fact dual cellular complexes. 
They are defined by a real positive function mg on the surface, the multiplicity. 

Consider a point x & S;as the set Q is discrete, the distance d{x, Q) is realized 
by geodesies of minimal length, generically a single one. Let mQ{x) e [1, oo) be 
the number of such geodesies. If mQ{x) = 1, there exists a vertex 7r(a;) S Q such 
that the shortest geodesic from x to ■k{x) is the only geodesic from x io Q with 
such a small length. 

The Voronoi 2-set associated to a vertex v in Q, is 7r~^(-y), that is to say the 
set of points of E closer to this vertex than to any other vertex in Q. Each 2-set 
of y is a connected component of m'^{\). 

Likewise, the 1-sets are the connected components of mQ^{2). They are asso- 
ciated to pairs of points in Q. 

The 0-sets are the connected components of mQ^([3, +oo)). Generically, they 
are associated to three points in Q. 

V IS a, cellular complex (see below) and the complex D is its dual (generically 
a triangulation), its vertices are the points in Q, its edges are segments {x,x') 
for x,x' £ Q such that there exist points equidistant and closer to them. 

Proposition 10. The Voronoi partition, on a closed Riemann surface with a 
flat metric with conic singularities, of a given discrete set of points Q containing 
the conic singularities, is a cellular complex. 
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Wc have to prove that 2-sets are homeomorphic to discs, 1-sets are segments 
and 0-sets are points. 

First, 2-sets are star-shaped, for every point x cfoser to v £ Q than to any 
other point in Q, along a unique portion of a geodesic, the whole segment [x, v] 
has the same property. 

2-sets are open, if x is closer to w G Q than to any other point in Q, as it is 
discrete, d(x, Q \v) — d{x, v) > 0. By triangular inequality, every point in the 
open ball of this radius centred at x is closer to v than to any other point in Q. 

So 2-sets are homeomorphic to discs. 

Let a; be a point in a 1-set. It is defined by exactly two portions of geodesies 
D,D' from x to y,y' € Q (they may coincide). By definition, the open sphere 
centred at x containing DUD' doesn't contain any point of Q so it can be lifted 
to the universal covering, where the usual rules of euclidean geometry tell us that 
the set of points equidistant to y and y' around a; is a submanifold of dimension 
1. As the surface is compact, if it is not a segment, it can only be a circle. Then, 
it's easy to see that the surface is homeomorphic to a 2-sphere and that y and 
y' are the only points in Q. But this is impossible because an euclidean metric 
on a 2-sphere has at least three conic singularities |Tro| . 

The same type of arguments shows that 0-sets are isolated points. □ 



Fig. 6: The Voronoi'/Delaunay decompositions associated to two points on a 
genus two surface. 



Proposition 11. Such Delaunay/Voronoi dual complexes are semi- critical maps 
of the surface. Hence any Riemann surface admits semi-critical maps. 

The edge in V dual to € Di is a segment of their mediatrix so is 

orthogonal to (x,a;'). Hence, equipped with the euclidean length on the edges, 
{y,D) is a semi-critical map.D 

Remark 5. Apart from Voronoi'/Delaunay maps, circle packings [CdV90| give 
another very large class of examples of interesting semi-critical decompositions 
(see Fig.IZl). 
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Fig. 7: Circle packing, the dual vertex to a face. 



The semi-criticality of a double map gives a coherent system of angles (f> in 
(0, tt) on the oriented edges of A. An edge {x, x') G Ai is the diagonal of a certain 
diamond; 4>{x,x') is the angle of that diamond at the vertex x. In particular, 
(p{x,x') ^ (p{x',x) a priori. They verify that for every diamond, the sum of the 
angles on the four directions of the two dual diagonals is 27r (see Fig. Then 
the conic angle at a vertex is given by the sum of the angles over the incident 
edges. 



y 



(x,x\ 



Fig. 8: A system of angles for a semi-critical map. 



3.2. Continuous limit. We state the main theorem, a converging sequence of 
discrete holomorphic functions on a refining sequence of semi-critical maps of 
the same Riemann surface, converges to a holomorphic function. Precisely: 

Theorem 3. Let S be a Riemann surface and {''A,£k)k£fi a sequence of semi- 
critical maps on it, with respect to the same flat metric with conic singularities. 
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Assume that the lengths Sk = S{''A) tend to zero and that the angles at the 
vertices of all the faces of the are in the interval [77, 27r — 77] with rj > 0. 

Let {fk)k&i be a sequence of discrete holomorphic functions fk G such 
that there exists a function f on S which verifies, for every converging sequence 
(xfe)feGM of points of S with each Xk G ''Aq, f{\imk{xk)) = limfc(/fc(xfe)), then 
the function f is holomorphic on S. 

Such a refining sequence is easy to produce (see Fig. [5]) but the theorem 
takes into account more general sequences. A more natural refining sequence, 
which mixes the two dual sequences is given by a series of tile centering proce- 
dures |GS87| : If one calls 0/2 the cellular decomposition constructed from by 
replacing each tile by four smaller ones of half its size, and r{(}/2), r*{(}/2) the 
double cellular decomposition it defines, one has 

r(0/2) = r(0)" U "r*(0) and the interesting following sequence: 

r{0) -^0^ r{0/2) ^ 0/2 ^ • • • ^ 0/2" ^ . . . 

/ \ / \ / \ . (3.1) 

r*(0) r*{<}/2) 

The horizontal arrows correspond to tile centering procedures, and the ascending, 
respectively descending arrows, to tile centering, resp. edge centering procedures. 
This sequence is not that exciting though since locally, the graph rapidly looks 
like a rectangular lattice. More interesting infiation rules staying at criticality 
can be considered too (see Fig. [2T|). 

The demonstration of the continuous limit theorem needs three lemmas: 




Fig. 9: Refining a semi-critical map. 



Lemma 2. Let {fk)k&'N be a sequence of functions on an open set Q C C such 
that there exists a function f on f2 verifying, for every converging sequence 
{xk)k£N of points of fl, f{[\m.k{xk)) = \\Tak{fk{xk)) ■ Then the function f is 
continuous and uniform limit of {fk) on any compact. 
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Taking a constant sequence of points, we see that (fk) converges to / point- 
wise. So with the notations of the theorem, {fk{xk)) converges to f{x) and 
{fi{xk))i(^n to f{xk). Combining the two, {f{xk)) converges to /(x) so / is con- 
tinuous. If the convergence was not uniform on a compac sett, then there would 
exist a converging sequence {xk) with {fk(xk) — f(xk)) not converging to zero. 
But / is continuous in x = \\m{xk) and {fk{xk)) converges to /(x), which, 
combined, contradicts the hypotheses. □ 

Lemma 3. Let (ABCD) be a four sided polygon of the Euclidean plane such 
that its diagonals are orthogonal and the vertices angles are in [ry, 27r — rj] with 
r] > 0. Let {M,M') be a pair of points on the polygon. There exists a path on 
(ABCD) from M to M' of minimal length £. Then 

MM' sin ri 
> -. 

£ - A 

It is a straightforward study of a several variables function. If the two points 
are on the same side, MM' = i and sinry < 1. If they are on adjacent sides, the 
extremal position with MM' fixed is when the triangle MM'P, with P the vertex 
of (ABCD) between them, is isocel. The angle in P being less than r?, > 
sin ^ > If the points are on opposite sides, the extremal configuration is 

given by Fig. M^-, where ^ = ^.D 




1. M, M' on adjacent sides. 2. M, M' on opposite sides. 

Fig. 10: The two extremal positions. 



Lemma 4. Let {A,£) be any double cellular decomposition and a G C^(<^) a 
closed l-form. The l-form f ■ a is closed for any holomorphic function f G f2(A) 
if and only if a is holomorphic. 

Just check. 

Proo/0 We interpolate each function fk from the discrete set of points ^Aq 
to a function fk of the whole surface, linearly on the edges of '^<^ and harmonicly 
in its faces. 

Let (Cfc) be a converging sequence of points in S. Each (k is in the adherence 
of a face of Let Xk,yk be the minimum and maximum of Re fk around the 
face. By the maximum principle for the harmonic function Re fk, 

Re fkixk) < Re fkiCk) < Re fkivk). 
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Moreover, the distance between Xk and C,k is at most 25k, as well as for y^. It 
implies that (xfc) and {yk) converge to x = lim(Cfc), {fk{xk)) and {fk{Vk)) to 
/(x), and (Re fkiCk)) to Re f{x); and similarly for its imaginary part. So, by 
Lemma [21 the function / is continuous, and is the uniform limit of (fk) on every 
compact set. In particular, it is bounded on any compact. 

By the theorem of inessential singularities, since / is continuous hence bounded 
on any compact set, and that conic singularities form a discrete set in S, to show 
that / is holomorphic, we can restrict ourselves to each element U C S oi a, eu- 
clidean atlas of the punctured surface (without conic singularities). We have an 
explicit coordinate z on [/. 

Let 7 be a homotopically trivial loop in U of finite length £. We are going to 
prove that fdz = 0. The theorem of Morera then states that / is holomorphic. 

Let us fix the integer k. By application of Lemma[3]on every face of crossed 
by 7, we construct a loop 7^ € Ci('^'C>), homotopic to 7, of length £{'jk) < 

(see Fig. [TT|) . As the diameter of a face of is at most 2Sk, all these faces are 
contained in the tubular neighbourhood of 7 of diameter ASk- Its area is ASki 
and it contains the set C of S between 7 and 7^^. 



Assume / is of class C-^, on the compact C, \df\ is bounded by a number M. 
Applying Stockes formula to fdz, 



So §^f{z)dz = \mi f{z)dz. Taking a sequence of class functions con- 
verging uniformly to / on C, we prove the same result for / simply continuous 
because all the paths into account are of bounded lengths. 




Fig. 11: The discretised path. 
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As (/fe) converges uniformly to / on C and the paths are of bounded lengths, 
we also have that {\ §^^{fk{z) — /(•z))rf-z|)fceN tends to zero. But because the 
interpolation is linear on edges of 'X), fk{z)dz = fkdZ, the second integral 

being the coupling between a 1-chain and a 1-cochain of But since fk and 
dZ are discrete holomorphic, fkdZ is a closed 1-form, and <f^^ fkdZ = 0. So 

.fk{z)dz tends to zero and 

f{z)dz = 0. 
□ 



3.3. Criticality. 

Proposition 12. Let a be a holomorphic l-form, f ■ a is holomorphic for any 
holomorphic function if and only if j^^ x)'^ ~ I{x' y') •^"'^ each pair of dual edges 
{x,x'), {y,y'). 

Let (x, y, x' , y') S <)2 be a face of <C>, the Cauchy-Riemann equation for / • a, 
on the couple {x,x') and {y,y') is the nullity of: 



iv.y') .{x,x'} 
-I- 



i{y,y') i{x,x') 



1 i!{x)-rf{v) f , f{x)+f{y') f , f{x') + f{y) f , f{x') + f{y') 



^(2/,y')' 2 7 2 J 2 J 2 

(y,x) (x.y') {y,x') (x',y') 

1 tf{x) + f{y) f , fix') + fiy) f , fix) + fiy') f , fix') + f{y') 



e{x,x')'^ 2 J 2 

(x,y) (y,x') ix,y') {y',x') 



(if.an) (y',x'} 

after having developed, used the holomorphy of a, then the holomorphy of /.□ 
So to be able to construct out of the holomorphic 1-forms dZ given by local 
flat isometrics, and a holomorphic function a holomorphic 1-form fdZ, we have 
to impose that for each face {x,y,x',y') G <)2, Z{x) — Z{y) ~ Z[y') — Z[x'). 
Geometrically, it means that each face of the graph is mapped by Z to a 
parallelogram in C. But as the diagonals of this parallelogram are orthogonal, it 
is a lozenge (or rhombus, or diamond). 

Definition 4. A double (yl, t) of a Riemann surface E is critical if it is semi- 
critical and each face of ()>2 are lozenges. Let 5{A) be the common length of their 
sides. 

Remark 6. This has an intrinsic meaning on Z", the faces of <^ are genuine 
lozenges on the surface and every edge of A can be realized by segments of 
length given by I, two dual edges being orthogonal segments. 
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Another equivalent way to look at criticality can be useful: a double {A, £) is 

critical if there exists an application Z : S \ P ^ C from the universal covering 
of the punctured surface E \ P for a finite set P C Aq into C identified to the 
oriented Euclidean plane such that 

— the image of an edge a G /li is a linear segment of length i{a), 

— two dual edges are mapped to a direct orthogonal basis, 

— Z is an embedding out of the vertices, 

— there exists a representation p of the fundamental group tti{S\P) into the 
group of isometrics of the plane respecting orientation such that, 

V7 e 7ri(i;\ p), Zo 7 = p{j) o z, 

— and the lengths of all the segments corresponding to the edges of <^ are all 
equal to the same 6 > 0. 

The criticality of a double map gives a coherent system of angles in (0, tt) 
on the unoriented edges of yl, (f>{x, x') is the angle in the lozenge for which {x, x') 
is a diagonal, at the vertex x (or x'). They verify that for every lozenge, the sum 
of the angles on the dual diagonals is n. Then the conic angle at a vertex is given 
by the sum of the angles over the incident edges. 

Every discrete conformal structure {A,£) defines a conformal structure on 
the associated topological surface by pasting lozenges together according to the 
combinatorial data (though most of the vertices will be conic singularities). Con- 
versely, 

Theorem 4. Every closed Riemann surface accepts a critical map. 

Proof^ We first produce critical maps for cylinders of any modulus: Consider 
a row of n squares and glue back its ends to obtain a cylinder, its modulus, the 
ratio of the square of the distance from top to bottom by its area is 

Stacking m such rows upon each other, one gets a cylinder of modulus ^. 

Squares can be bent into lozenges yielding a continuous family of cylinders of 
moduli ranging from zero to ^ (see Fig. [T^ . Hence we can get cylinders of any 
modulus. 



\\\\\\\-\ 



Fig. 12: Two bent rows. 



Dehn twists can be performed on these critical cylinders, see Fig. [131 
Gluing three cylinders together along their bottom [n has to be even), one 
can produce trinions of any modulus (see Fig. I14p and these trinions can be 
glued together according to any angle. Hence, every Riemann surface can be so 
produced [Bus] . □ 
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Fig. 13: Performing a Dchn twist. 




Fig. 14: Gluing three cylinders into a trinion. 

Remark 7. An equilateral surface is a Riemann surface which can be triangulated 
by equilateral triangles with respect to a flat metric with conic singularities. 
Equilateral surfaces are the algebraic curves over Q |VoSh| so are dense among 
the Riemann surfaces. Cutting every equilateral triangle into nine, three times 
smaller, triangles (see Fig. [TS)) . one can couple these triangles by pairs so that 
they form lozenges, hence a critical map. 

In Figures I16H19I are some examples of critical decompositions of the plane. 
In Fig. [20l a higher genus example, found in Coxeter [Coxl] . of the cellular 
decomposition of a collection of handlebodies (the genus depends on how the 
sides are glued pairwise) by ten regular pentagons, the centre is a branched 
point of order three; together with its dual, they form a critical map. It is the 
case for any cellular decomposition by just one regular tile when its vertices are 
co-cyclic. This decomposition gives rise to a critical sequence using the Penrose 
inflation rule |GS87j . Fig. [2T] illustrates this inflation rule sequence on a simpler 
genus two example where each outer side has to be glued with the other parallel 
side. 



3.4- Physical interpretation. 

Theorem 5. A translationally invariant discrete conformal structure (yl, p) on 
A the double square or triangular /hexagonal lattices decomposition of the plane 
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Fig. 15: An equilateral triangle cut in nine yielding lozenges. 




Fig. 16: A 1-parameter family of critical deformations of the square lattice. 

or the genus one torus, is critical and flat if and only if the Ising model defined 
by the interaction constants '■= iArcsinhpe on each edge e & Ai is critical as 
usually defined in statistical mechanics \McCWi . 

Proof . We prove it by solving another problem which contains these two par- 
ticular cases, namely the translationally invariant square lattice with period 
two |Yamj . At a particular vertex, the flat critical condition on the four confor- 
mal parameters is: 

4 

arctan pi = tt, 

1=1 



which is obviously invariant by all the symmetries of the problem, including 
duality. When pi ~ Pi+2, we get the usual period one Ising model criticality on 
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Fig. 17: A 2-parameters family of critical deformations of the triangu- 
lar/hexagonal lattices. This family, key to the solution of the triangular Ising 
model, induced Baxter to set up the Yang-Baxter equation |Bax| . Our notion of 
criticality fits beautifully into this framework. 




Fig. 18: The order 5 Penrose quasi crystal. 

the square lattice 

sinh2i^/iSinh2/^„ = 1, 

and likewise when one of the four parameters degenerates to zero or infinity, the 
three remaining coefficients fulfill 

sinh 2Ki sinh 2Kji sinh 2Kiii = sinh 2Ki + sinh 2Kii + sinh 2Kiii 

which is (a form of) the criticality condition for the triangular/hexagonal Ising 
model. The case shown in Fig. [TB] occurs when pi = pj, = 1, implying P2Pa = 1- 
□ 
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Fig. 19: Lozenge patchworks. 




Fig. 20: Higher genus critical handlebody. 



We see here that flat criticality, when the angles at conic singularitites are 
multiples of 27r, is more meaningful than criticality in general. This theorem 
is important because it shows that statistical criticality is meaningful even at 
the finite size level. It is well known [KW| that for lattices, it corresponds to 
self-duality, which has a meaning for finite systems; here we see that self-duality 
corresponds to a compatibility with holomorphy. In a sense, our notion of crit- 
icality defines self-duality for more complex graphs than lattices. Furthermore, 
we will see in Sect. |4] that criticality implies the existence of a discrete mass- 
less Dirac spinor, which is the core of the Ising model. Although wc saw that 
criticality implies a continuous limit theorem, the thermodynamic limit is not 
necessary for criticality to be detected, and to have an interesting meaning. 

It is easy to produce higher genus flat critical maps and compute their crit- 
ical temperature, the examples in Figures [20ll2T] have four kinds of interactions 
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Fig. 21: Sequence of critical maps of a genus two handlebody using Penrose 
inflation rule. 



corresponding to the diagonals of the two kinds of quadrilateral tiles. They are 
critical when the angles of the quadrilaterals are j,^,^, and corresponding 
to Ising interactions 

TL7T 

sinh 2 A'„ = tan — . (3.2) 
10 ^ ' 

The author had made no attempt to verify these values numerically. 

A general way is, considering a critical genus one torus made up of a transla- 
tionally invariant lattice, to cut two parallel segments of equal length and seam 
them back, interchanging their sides. This creates two conic singularities where 
an extra curvature of —2tt is concentrated at each point, yielding a genus two 
handlebody. Repeating the process, we may produce critical handlcbodics of ar- 
bitrarily large genus if we start with a very fine mesh. One has to beware that 
our continuous limit theorem applies only to fixed genus, it cannot grow with 
the refinement of the mesh. This explains why the union-jack lattice (the square 
lattice and its diagonals) or the three dimensional Ising model, which can be 
modelled as a genus mnp surface for a 2m x 2n x 2p cubic network, are beyond 
the scope of our technique as far as a continuous limit theorem is concerned. 
With this restriction in mind, we see that both the existence and the value of 
a critical temperature is essentially a local property and neither depends on the 
genus nor on the modulus of the handlebody. It is not the case for more interest- 
ing quantities such as the partition function, which can be obtained in principle 
from the discrete Dirac spinor that criticality provides, defined in Sect. S) But 
such a calculus is beyond the scope of this article. 

Apart from the standard lattices, the critical temperature of other well known 
graphs can be computed using our method, for example the labyrinth |BGB| . 
whose diamond is pictured in Fig. [22l has the topology of the square lattice but 
has five different interactions strengths controlled by two binary words, labelling 
the columns and rows by O's and I's. And also new ones such as the "street 
graph" depicted in Fig. [23l Its double row transfer matrix appears to be the 
product of three commuting transfer matrices, two triangular and a square one. 

Other cases such as the Kagome |Syo| or more generally lattices of chequered 
type |Utij can be handled using a technique called electrical moves |CdV96j 
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22: The diamond graph of a critical labyrinth lattice. 




Fig. 23: The "street" lattice. 
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which enables us to move around, and causes appearing or disappearing conic 
singularities of a flat metric. This will be the subject of a subsequent article, 
explaining the relationship between discrete holomorphy, electrical moves and 
knots and links. These electrical moves act in the space of all the graphs with 
discrete conformal structures in a similar way to that of the Baxterisation pro- 
cesses in the spectral parameter space of an integrable model (see [AdABM] ). 
We are going to see that the link with statistical mechanics is even deeper than 
simply pointing out a submanifold of critical systems inside the huge space of all 
Ising models, as the similarity with the continuous case extends to the existence 
of a discrete Dirac spinor near criticality. 

3.5. Polynomial ring. 

Definition 5. Let {A,£) be a critical map. In a given flat map Z : U C on the 

simply connected U, choose a vertex zq G Aq, and for a holomorphic function f , 
define the holomorphic functions and f by the following formulae: 

f\z) ■.^e{z)f{z), 

where f denotes the complex conjugate and s{r) = = —1, 

/'W:4(£/..Z)\ 

See |Duf| for similar definitions. Notice that /' is defined up to e if one changes 
the base point. 

Proposition 13. Let {A,i)be a critical map. In a given flat map Z : U ^ C on 
the simply connected U, for every holomorphic function f e df ~ f'dZ. 

We hence call f the derivative of f . 

Consider an edge {x, y) S <^i, x G Iq, y <E F^, 

So = -^^^^i2(s) - Z{xmZ(y) - Z(x)) = /(a) - /(!).□ 

Definition 6. Let U be a simply connected flat region and zq £ U . Define in- 
ductively the holomorphic functions Z^{z) := J^^ ^Z^~^dZ given Z^ := 1. As 
the space of holomorphic functions on U is finite dimensional, these functions 
are not free; let Pu be the minimal polynomial such that Pu{Z) = Z" + . . . = 0. 

Conjecture 6 The space of holomorphic functions on U , convex, is isomorphic 
to C[Z]/Pu. 
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Wc won't define here the notion of convexity, see [CdV96| . The question is 
whether the set {Z'') generates the space of holomorphic functions. The problem 
is that zeros are not locahsed, and as the power of Z'^ increases, the set of its 
zeros spread on the plane and get out of U. Figure [24] is an example on the unit 
square lattice with U the square [—10, 10] © [—10, 10]i, the degree increases with 
k until 16 where four zeros get out of the square. So a definition of the degree 
of a function by a Gauss formula is delicate. 



Z^^ and its zeros. Z^^ and its zeros. 

Fig. 24: The zeros of Z'^^ get out of the square [-10, 10] ® [-10, lQ]i. 

4. Dirac Equation 

Although we believe our theory can be applied to a lot of different problems, our 
motivation was to shed new light on statistical mechanics and the Ising model 
in particular. This statistical model has been linked with Dirac spinors since the 



work of Kaufman 



K] and Onsager and Kaufman [KOj . We refer among others 



to [McCW81|ISMJ|IKC] . Hence we are interested in setting up a Dirac equation in 



the context of discrete holomorphy. To achieve this goal we first have to define 
the discrete analogue of the fibre bundle on which spinors live. We therefore 
have to define a discrete spin structure. Physics provides us with a geometric 
definition [KC| based on paths in a certain Z2-homology, that we generalise to 
our need (higher genus, boundary, arbitrary topology). We begin by showing 
that such an object in the continuum is indeed a spin structure, then define the 
discrete object. We then set up the Dirac equation for discrete spinors, show 
that it implies holomorphy and that the existence of a solution is equivalent to 
criticality. The Ising model gives us an object which satisfies the discrete Dirac 
equation, namely the fermion, W = cr/i as defined in [KG] , corresponding to a 
similar object defined previously by Kaufman [K]. It fulfills the Dirac equation 
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at criticality, but also off criticality, corresponding to a massive Dirac spinor. 
We will end this article by describing ofF-criticality, as defined by the author's 
Ph.D. advisor, Daniel Bennequin. 

4.1. Universal spin structure. A spin structure [Mil] on a principal fibre bundle 
(E, B) over a manifold B, with SO(n) as a structural group, is a principal fibre 
bundle {E',B), of structural group Spin(n), and a map f : E' ^ E such that 
the following diagram is commutative: 

E' X Spin(7i) -> E' 

\ 

i /X A if B 

/ 

E X S0(7i) E 

where A is the standard 2-fold covering homomorphism from Spin(7T,) to SO(n). 

In this paper we consider only spin structures on the tangent bundle of a 
surface. On a generic Riemann surface 17, there is not a canonical spin structure. 
We are going to describe a surface 17, 2^~'=^(^)-fold covering of 17, on which there 
exists a preferred spin structure. It allows us to define every spin structure on 17 
as a quotient of this universal spin structure. We will treat the continuous case 
and then the discrete case. 

Definition 7. Let S be a dijjerentiable surface with a base-point y"; E is the 
set of pairs (z, [A]2), where z £ S is a point and [A]2 the homology of a path A 
from to z considered in the relative homology Hi{E, {y'^, z}) (Xi Z2. 

E is the 2^"^'^' covering associated to the intersection H of the kernels of 
all the homomorphisms from 7ri(I7) to Z2, that is to say the quotient of the 
universal covering by the subgroup H C t^i{E) of loops whose homology is null 
modulo two. 

Choose vq a tangent vector at For each point z € E, define E^ := E \ 
{y^, z} U U the blown up of E at j/° and z (add only one circle in the case 
y*^ — z). Consider the set of oriented paths in E^, from the point corresponding 
to the vector vq at to the directions at z (the vector vg is needed only when 
z = Define an equivalence relation ~2 (see Fig. [25| on this set by stating 
that two paths A, A' are equivalent if and only if A — A' is a cycle and [A — A'] 2 = 
in the homology Hi{E \ {z}, Z2). 

Definition 8. The universal spin structure S of E is the set of pairs (z, [A]^^), 
with z ^ E and [A]^, the ^^-equivalence class of the path A from ?/° to z in E^- 

Theorem 7. S is a spin structure on E and is the only one such that the action 
of the fundamental group t^i^E) on E can be lifted to. Moreover it is the pull-back 
of any spin structure on E . 

Proof . The proof is in three steps, we check that 5 is a spin structure, we define 
a spin structure Sq through group theory and we show that both are equal to a 
third spin structure Si. 
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Fig. 25: Paths of different classes with respect to for z ^ -ip and z = y". 



There is an obvious projection from S to 12 defined by (z, [A]^^) ^ (z, [A]2). 
The fibre of this projection at (z, [A]2) is the set of ^^-equivalence classes of 
paths from y° to the blown-up circle at z. To each class is associated the tangent 
direction at z so is a covering of ST,XS). As Hi{I] \ {z},Z2) is 23-x(-S') 
dimensional (a loop around z is not homologically trivial), for each point in 
ST{IJ), there are two different lifts. The path in Sz corresponding to turning 
around z once yields the Z2-deck transformation. Hence 5 is a spin structure on 

Let G := 7ri(Z') and G' := tti{STS); the §i-fibre bundle ST{S) S induces 
a short exact sequence 1, ^ G' ^G. Every double covering of ST 17 is defined 
by the kernel S" of an homomorphism u from G' to Z/2, moreover, for 5" to be 
a spin structure, its intersection with the subgroup Z must be 2Z. 

Likewise, the fibration S ^ S implies that the fundamental group H' := 
TTi{STS) of the directions bundle of IJ is the subgroup of G" over H := 



Z - 


^ H' - 


^ H 




i 


i 


i , 


(4.1) 


Z - 


^ G' - 


G 





The intersection of the subgroups H' and S' is a well defined spin structure 
5o on S: Indeed, consider another spin structure S" = Ker [v : G' ^ ^/2) on 
17, its intersection with Z is 2Z hence the kernel oi u — v contains the whole 
subgroup Z, that is to say u — v comes from a homomorphism of G to Z/2 and 
we have S" H H' = S' C] H' . In other words, Sq is the unique spin structure on 
S which is the pull-back of a spin structure on S and it is the pull-back of any 
spin structure. 

Let z € 17 be a point, consider the set of paths in STS from the base point 
(y'',v°) to any direction at z. Consider on this set the equivalence relation 
defined by fixed extremities Z/2-homology. The class [A]^' of a path A from 
{y°,v°) to {z,v) is its homology class in Hi{STE , {{y° ,v°)^{z,v)}) ® Z/2. The 
projection STU splits Hi{STE, {(y°, v°), {z,v)}) (g) Z/2 into 

Z/2 ^ HiiSTS, {(2;°,t^°), (z, v)}) ® Z/2 ^ Hi{S, {y", z}) ® Z/2, (4.2) 

hence the set Si of pairs (z, [A]^' ) for all points z lE S and all paths A, is a spin 
structure on E. 
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Let 5" be a spin structure on S, it defines an element in Z/2 for eaclr loop 
in STS . So each path in STS beginning at (y'^, v^) defines, through the split- 
ting [421 an element in Si which is then the pull-back of S" to ^ , hence Sq ~ Si. 

On the other hand S = Si because there is a continuous projection from S to 
Si: For an element (z, [A]^,), consider a C^-path A S 17 representing the class. 
Lift it to a path in STS by the tangent direction at each point, its class [A]^'_ 
only depends on [A]^_ and gives us an element in 5i. □ 



4-. 2. Discrete spin structure. 

Definition 9. Let Y be a cellular complex of dimension two, a spin structure 
on T is a graph T' , double cover of the l-skeleton of T such that the lift of the 
boundary of every face is a non-trivial double cover. They are considered up to 
isomorphisms. Let So be the set of such spin structures. 

A spinor tjj on T' is an equivariant complex function on T' regarding the 
action o/ Z/2, that is to say, for all ^ G Tq, = — */ C represents the 
other lift. 

Remark 8. Usually, a spinor field is a section of a spinor bundle, that is to say a 
square root of a tangent vector field. Here, we consider square roots of covectors; 
we should say cospinors. 

A discrete spin structure is encoded by a representation of the cycles of T, 
Zi{T) := Kcr d n Ci(T), into Z/2 which associates to 7 G Zi{T), the value 
m(7) = if it can be lifted in Y' to a cycle and /i(7) = 1 if it can not. By 
construction, the value of the boundary of a face is 1 and the value of a cycle 
which is the boundary of a 2-chain of T is the number of faces enclosed, modulo 
two. 

We are going to show that this structure is indeed a good notion of discrete 
spin structure. First, there are as many discrete spin structures on a surface as 
there are in the continuous case: 

Proposition 14. On a closed connected oriented genus g surface S , the set Sd 
of inequivalent discrete spin structures of a cellular decomposition T is of cardinal 
2^^. The space of representations of the fundamental group of the surface into 
Z/2 acts freely and transitively on Sd- 

We explicitly build discrete spin structures and count them: Let T be a maximal 
tree of T, that is to say a sub-complex of dimension one containing all the vertices 
of T and a maximal subset of its edges such that there is no cycle in T. Choose 
2g edges {ek)i<k<2g in T \ T such that the 2g cycles (7^) G Zi{T)'^3 extracted 
from (T U ek)i<k<2g form a basis of the fundamental group of S (and T). Let 
T+ := TUk ek and consider T', the sub-complex of the dual T* formed by all the 
edges in Y* not crossed by T+. It is a maximal tree of Y* . Likewise we define 
n:=rUfee*. , , , 

We construct inductively a spin structure Y': its first elements are a double 
copy of T and we add edges without any choice to make as we take leaves out of 
T'j^. When only cycles are left, a choice concerning an edge has to be taken, 
opening a cycle in T'j^. The process goes on until T'j^ is empty. 
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These choices are completely encoded by a representation fi such as in the 
remark, and the 2g values {fJ-{'jk))i<k<2g determine the spin structure. On the 
other hand, this representation defines the spin structure and there are 2^^ such 
different representations. Hence the choices of the maximal tree and the edges 
Bk are irrelevant. 

Because a cycle in T belongs to a class in the fundamental group of the 
surface (up to a choice of a path to the base point, irrelevant for our matter), 
the representations of the fundamental group into Z/2 obviously act on spin 
structures: A representation p : 7Ti{S) Z/2 associates to a spin structure 
defined by a representation fi : Zi{T) — > Z/2, the spin structure defined by the 
representation p{^) such that p{p){j) := /i(7) + ^([7]), where [7] G 7i'i(^) is the 
class of the cycle 7 in the fundamental group. This action is clearly free, and 
transitive because the set of representations is of cardinal 2^^. □ 

Given A ~ F Li F* a double cellular decomposition, we introduce a cellular 
decomposition which is the discretised version of the tangent directions bundle 
of both F and F*: 

Definition 10. The triple graph T is a cellular complex whose vertices are un- 
oriented edges of(), Tq = {{x , y\ / (x , y) e ^1}. Two vertices {x,y\,{x' ,y'} G To 
are neighbours in T iff the edges {x,y) and {x',y') are incident (that is to say 
X = x' or X — y' or y = x' or y — y' ), and they hound a common face of <0>. 
There are two edges in T for each edge in A. For this to he a cellular decompo- 
sition of the surface in the empty boundary case, one needs to add faces of three 
types, centred on vertices of F , of F* and on faces of (see Fig. \26\} . 




Fig. 26: The triple graph T. 



Remark 9. The topology of the usual tangent directions bundle is not at all 
mimicked by the incidence relations of T, the former is 3 dimensional and the 
latter is a 2-cellular complex. 

Let (a;°, y°) £ be a given edge. All the complexes F, r"*, <>, T are lifted to 
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Definition 11. The discrete universal spin structure T' is the following 
1-complex: Its vertices are of the form {{x,y), [7^0]), where {x,y) ^Tq is a pair 
of neighbours in <C> and 7^0 is a path from to y on F* , avoiding the faces x* 
and x'^* . We are interested only in its relative homology class modulo two, that 
is to say [7^^] G \ a;*, {y", y}) (8) Z2. We will denote a point by ((a;, y), 7^0) 

and identify it with ((x, y), 7'^o) whenever 7^0 and 7'^o o-fs homologous. 
Two points ((a;,J/),7^o) and ((a;', ?/'), 7^, ) are neighbours in T' if 

— x = x' , {y, y') S and j-j" — 7^, + (y, y') is homologous to zero in Hi{r* \ 
a;*)(8)Z2, 

— y — y' , (x, x') <E Ti and 7^ — 7^, is homologous to zero in Hi (F* \ x*) ® Z2 . 

T' is a double covering of T and it is connected around each face (see Fig.[27|). 
It is a discrete spin structure on T in the sense defined above. Once a basis of the 
fundamental group tti (T) is chosen, every representation of the homology group 
of U into Z2 allows us to quotient this universal spin structure into a double 
covering of T, yielding a usual spin structure T' . 



4-. 3. Dirac equation. A spinor changes sign between the two lifts in T' of a 
vertex of T, in other words it is multiplied by —1 when it turns around a 
face. The faces of Y which are centred on diamonds arc four sided. We set 
up the spin symmetry equation for a function C, on Tq, on a positively ori- 
ented face (Ci, ^2,^31 ^4) £ ^2 around a diamond, lifted to an 8-term cycle 




Fig. 27: Double covering around faces of T. 




c(e.+). 



(4.3) 
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The coherent system of angles (p given by a semi-critical structure locally 
provides a spinor respecting the spin symmetry away from conic singularities: 
Define half angles 9 on oriented edges of T in the following way: Each edge 
(C, £.') e Ti cuts an edge a € Ai, set 6'(^, := ±^ whether turns in the 

positive or negative direction around the diamond. Choose a base point G Tq, 
define C, by C('^o) = 1 and 

C(0 :=expz^0(A) (4.4) 

for any path 7 from ^0 to ^. The sum of the half angles are equal to tt around the 
faces of <^ and half the conic angle around a vertex, so if it is a regular flat point, 
we get ^ = TT again, hence C is a well defined spinor. As diagonals of the faces of 
<^ are orthogonal, C fulfills the spin symmetry. Moreover, if the conic angles are 
congruous to 27r modulo 47r, C can be extended to any simply connected region; 
if the fundamental group acts by translations, C is defined on the whole T' . 

We are going to define a propagation equation which comes from the Ising 
model. It is fulfilled by the fermion defined by Kaufman [K] which is known to 
converge to a Dirac spinor near criticality. We will use the definition ip = crfi 
given by Kadanoff and Ceva [KC| . The Dirac equation has a long history in 
the Ising model, beginning with the work of Kaufman [IT and Onsager and 
Kaufman |KOj . we refer among others to [McCW Sl , SMJ , KCj . The equation 
that we need is defined explicitly in |DD| , hence we will name it the Dotsenko 
equation, even though it might be found elsewhere in other forms. It is fulfilled 
by the fermion at criticality as well as off criticality. But this equation is only a 
part of the full Dirac equation. For a function ^ on Tq, with the same notations 
as before, and if a e /li is the diagonal of the diamond, between (C2,C3) and. 
(C4,Ci) (see Fig. EH): 

a^t) = V^T^a^t) - pia)aci)- (4.5) 

A check around the diamond shows that it also implies that C is a spinor: We 
write the Dotsenko equation in and , 

cm =v/ l + p(a*m 3+) - Piaiaa), 
hence, as y^l + p(a)2-^l + p{a*Y = p(a) + p{a*), 

ait) ^p{a*)a£,t)-v^+^m) 
= -c(er)- 

The Dirac equation is the conjunction of the symmetry (|4.3[) and the Dot- 
senko (|4.5p equations. We will see that this same equation describes the massive 
and massless Dirac equation, the mass measuring the distance from criticality. 

Given two spinors C,' , their pointwise product is no longer a spinor but a 
regular function on T . As there are two edges in T for each edge in A, there 
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Fig. 28: The Dotsenko equation. 



is an obvious averaging map from 1-forms on Y to 1-forms on A: Wc define 
drCC € C'^(yl) by the foUowing formula, with the same notation as before, 

2 / drcc c(6)c'(6) - c(6)c'(6) + auK'iCi) ~ c(a)c'(Ci)- 

J a 

dxCC is by definition an exact 1-form on Y but its average is not a priori exaet 
on A. 




Fig. 29: The 1-form on A associated to two spinors. 



Proposition 15. If and C respect whether the spin symmetry or the Dotsenko 
equation, then drCC a closed 1-form. If C, is a Dirac spinor and Q' fulfills the 
Dotsenko equation, then drCC holomorphic, dxCC anti-holomorphic and every 
holomorphic 1-form on A can be written this way on a simply connected domain, 
uniquely up to a constant. 

A sufficient condition for drCC to be closed on A is that, with the same 
notations as above, C(e3)C'(6) " C(6)C'(6) = C(^4)C'(^4) - C(6)C'(ei) because 
iay* drCC for a vertex y £ Aq is a sum of such differences on the edges of Y 
around y. This is so if there exists a 2 x 2-matrix A such that 

Wt)) - ^ Wt)) ' 
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a similar formula for C and *^ ( q '^j. ) ^ ~ (o ^1 J ' "^^^ solutions are of 



the form A = [ ^ ^ ^ ^ , — -— ttt I for a complex number A G C, e = ±1 and a 



determination of v 1 + A^. This is the case for the spin symmetry, A = — i,e = +l 
and for the Dotsenko equation, A = p{a), e = — 1, vT+A^ > 0. 
If C is a Dirac spinor and C fulfills the Dotsenko equation, then 



=<fe+)( v/i+pH Vfe+) - p(«)c'(e2+)) 
=*p(«)(c(e3^)c'(C3+)-c(C2+)c'(e2+)) / drCC. 

J a 

So drCC is holomorphic. Of course, d((' is anti-holomorphic. Conversely, if drCC 
is holomorphic with C a Dirac spinor, then ^' fulfills the Dotsenko equation. 

Given a holomorphic 1-form a G i?'"'^'^-' (yl), define ar on Ti by the obvious 
map f,, , f ar '■= L ,\ a. It is a closed 1-form on T because a is closed 
on A, so there exists a function a on any simply connected domain of Tq, unique 
up to an additive constant, such that dra = ar- A check shows that the only 
spinors Q" such that drCC" = OonA are the one proportional to C. It is consistent 
with the fact that the Dirac spinor is of constant modulus (see Eq. (|4.6p ). Hence 
the function C := a/C on T' is the unique spinor (up to a constant times 1/C ~ C) 
such that drCC ~ ct-O 

Notice that for ( a Dirac spinor, the holomorphic 1-form associated to it on 
A is locally, for a given flat coordinate Z, drCC = XdZ, with A G C a certain 
constant. 



4-.4- Existence of a Dirac spinor. 

Theorem 8. There exists a Dirac spinor on a double map iff it is critical for 
a given flat metric with conic angles congruous to 2tt modulo Att and such that 
the fundamental group acts by translations. The Dirac spinor is unique up to a 
multiplicative constant. 

Proof\^ Let be a non-zero Dirac spinor. Consider a positively oriented face 
(Cij'C2,'C3,C4) G ^2 around a diamond with diagonals a, a* as in Fig. [551 lifted to 
an 8-term cycle 

, ^2+ , > > ^1 . ^2 > ^3 ,^i)eC,{r). The equation 

^i±<£l ^ p{a*)+i 
defines an angle 0(a) G (0, tt) for every edge a d Ai. 
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The Dotscnko and symmetry equations combine into 

C(e^) = -7=rhci^t)- (4-6) 

The fact that ^ is a spinor impHes that, summing the four angles around the dia- 
mond, we get e'('A(a)+0(o )) = _i ^s each angle is less than tt, their sum is equal 

to TT. The same consideration around a vertex x £ Ao, yields exp i J2{x x')eAi '^^^^ ' 
— 1. So is a coherent system of angles and the map is critical with conic angles 
congruous to 27r modulo 47r. 

Conversely, given a coherent system of angles with conic angles congruous 
to 27r modulo 47r, the preceding construction described by Eq. (|4.4p gives the 
only Dirac spinor. □ 

In this case, dZ is a well defined holomorphic 1-forni on the whole surface. 

Corollary 9 Let (yl, p) he a discrete conformal structure and P a set of ver- 
tices, containing among others the vertices v such that the sum X]e~ii Arctanp(e), 
summed over all edges e incident to v, is greater than 2tt. The discrete conformal 
structure is critical with P as conic singularities if and only if there exist Dirac 
spinors on every simply connected domain containing no point of P. 

We define in which sense a discrete spinor converges to a continuous spinor. 
We don't define these spinors on specific spin structures but rather on the uni- 
versal spin structure S. 

Consider a sequence of finer and finer critical maps such as in Theorem |3l 
Choose a converging sequence of base points j/^) G '^Tq on each critical map 

such that the direction sequence {-j^r^) converges to a tangent vector (a;°, v^). 

Consider a sequence of points {xk,yk) G '"^O: defining a sequence of points 
(xk) converging to x in 17 and a converging sequence of directions v = lim^^|^^^^^. 
By compacity of the circle, there exist such sequences for every point x £ S 
and the criticality implies that it is in at least three directions for fiat points, 
separated by angles less than tt. 

The different limits allow us to identify, after a certain rank, the relative 
homology groups i/i \ x* , {y^ , } ) ® Z2 with ili (r„ { (a;° , ) , (x, i;) } ) ® Za , 
the classes of paths in the blown-up of S at x° and x. 

Definition 12. We will say that a sequence {Ck)ki£N of spinors converges if and 
only if, for any converging sequences, {{xk,yk) S '^^0) k^f^ defining a limit tan- 
gent vector, and {[Xk])k£n of classes of paths in 'T* from to y, avoiding the 
face x^, the sequence of values {Ckixk, [^k])) converges. 

Remark 10. It defines a continuous limit spinor C, by equivariance: Let x £ S, 
the set Dx of directions in which there exist converging sequences of discrete 
directions is by definition a closed set. Let u, v two boundary directions of 
such that the entire arc A of directions between them is not in D^. Consider 
[(x, [X]x), (x, [A']a;)] C S a. lift of A. The circle §^ acts on the directions, hence on 
the ~a;-classes, let ■0 e (0, tt) the angle such that (x, e*'''[A]a;) = (x, [A'J^). Define 

C(x,e^^[A],):=e™(^)C(x, [A],), 
where = (^(</') — f for Dirac spinors). 
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Theorem 10. Given a sequence of critical maps such as in Theorem with 
Dirac spinors on all of them, they can be normed so that they converge to the 
usual Dirac spinor on the Riemann surface. 

In a local flat map Z, the square of the discrete Dirac spinor on '^T' is (up to 
a multiplicative constant) the 1-form dZ evaluated on the edges. Hence their 
sequence converges. 



4-.5. Massive Dirac equation, discrete fusion algebra and conclusions. For com- 
pleteness and motivation, wc describe below the situation off-criticality where 
elliptic integrals come into play, and investigate a form of the discrete fusion 
algebra in the Ising model. This work was done by Daniel Bennequin and will 
be the subject of a subsequent article. 

A massive system in the continuous theory is no longer conformal. In the 
same fashion, Daniel Bennequin defined a massive discrete system of modulus k 
as a discrete double graph {A, p) such that, for each pair (a, a*) of dual edges, 

p{a)p{a*)^^-. (4.7) 

The massless case corresponds to A: = 1. We showed that criticality was 
equivalent to a coherent system of angles (j){a) such as shown in Fig. [8l defined 
by tan = p(a), and adding up to 27r at each vertex of the double, except at 

conic singularities. The Dirac spinor was constructed using the half angles 
Similarly, for every edge, we define the massive "half angle" uia) as the elliptic 
integral 

where the measure is deformed by 

k^ + k''^ = 1, (4.9) 



Z\((^) := a/I - Psin^, (4.10) 



A'{if) := ^1-k'^ If. (4.11) 

Using these non-circular half angles, and the corresponding "square angle" 
Ik' /q^ A'^ip) ' ^"^^ "^^'^ construct a massive Dirac spinor wherever the following 
"flatness" condition is fulflUcd: 



{Ik' - u{a)) = Ik' mod Alk' for each face F <G A2, (4.12) 
^^(u(a)) = Ik' mod Alk' for each vertex v S Aq. (4.13) 



aedF 



Daniel Bennequin noticed that the fusion algebra of the Ising model could be 
understood at the flnite level: Consider a trinion made of cylinders of a square 
lattice, of width m and n, glued into a cylinder of width to+ti. It has been known 
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since Kaufman [K] that, in the transfer matrix description of the Ising model, 
the configuration space of the Ising model on each of the three boundaries is a 
representation of spin groups spin(?Ti), spin(n) and spin(m + 7i) respectively. If m 
is odd, there exists a unique irreducible representation A of spin(m) but when 
m is even, there are two irreducible representations, Zi+ and A~ . A pair of pants 
gives us a map spin(m) x spin(ri) spin(m + n), in the case of a pair of pants of 
height zero, it's the inclusion given by the usual product. The representations of 
spin(m + n) induce representations of the product group that can be split into 
irreducible representations. If the three numbers are even, 

A+ A+ (E) A+ + A- A- , (4.14) 
A- ^ A+ ® A- + A- (E) A+, (4.15) 

(4.16) 

while if only one of them is even, 

A^ A+ (EA + A- (EA, (4.17) 

(4.18) 

A+~^A(E)A, (4.19) 
A-^A®A. (4.20) 

Let us compile these data in an array and relabel A by a, by 1 and A^ by 



and if m and n are both odd. 





1 


e 


a 


1 


1 


e 


a 


e 


e 


1 


a 


a 


a 


a 


1 + e 



(4.21) 



This is read as follows, the 1 + e in the slot a ® a ioj example, means that the 
representation 1 and the representation e of spin(m + n) both induce a factor 
a ® a va. the representation in the product group spin(m) x spin(n). 

We get exactly the fusion rules of the Ising model. The only difference com- 
pared with the continuous case is that the algebra is not closed at a finite level. 
The columns, rows and entries are not representations of the same group, rather 
we have a product of representations of spin(n) and spin(m) as a factor of a 
representation of spin(n + m). 

These results provide evidence that a discrete conformal field theory might 
be looked for: the discrete Dirac spinor at criticality is the discrete version of 
the conformal block associated with the field ^ and some sort of fusion algebra 
can be identified at the finite level. The program we contemplate is, first to 
investigate other statistical models and see if there are such patterns. If that is 
the case, we must then mimic in the discrete setup the vertex operator algebra 
of the continuous conformal theory. This can be attempted by defining a discrete 
operator algebra, in a similar fashion to Kadanoff and Ceva [KC| , and splitting 
this algebra according to its discrete holomorphic and anti-holomorphic parts. 
The hope is that some aspects of the powerful results and techniques defined by 
Belavin, Polyakov and Zamolodchikov [BPZ| will still hold. A very interesting 
issue would be, as we have done for the Ising model, to realize the fusion rules 
of a theory in the discrete setup, yielding its Verlinde algebra. 
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